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When a dense stream of dry, non-cohesive grains hits a Hxed target, a collimated sheet is ejected 
from the impact region, very similar to what happens for a stream of water. In this study, as a 
continuation of the investigation why such remarkably different incident fluids produce such similar 
ejecta, we use discrete particle simulations to collide two unequal-width granular jets head-on in 
two dimensions. In addition to the familiar coherent ejecta, we observe that the impact produces 
a far less familiar quasi-steady-state corresponding to a uniformly translating free surface and flow 
held. Upon repeating such impacts with multiple continuum huid simulations, we show that this 
translational speed is controlled only by the total energy dissipation rate to the power 1.5, and is 
independent of the details of the jet composition. Our hndings, together with those from impacts 
against hxed targets, challenge the principle of scattering in which material composition is inferred 
from observing the ejecta produced during impact. 


I. INTRODUCTION 

If we collide two ideal fluid streams head-on in two di¬ 
mensions, one larger than the other but both the same 
speed, we would observe a steady state with two equally 
sized and inclined streams leaving the impact region. 
This is because the larger jet has more incident mo¬ 
mentum than the smaller, and these ejecta are how the 
system accommodates for the momentum excess. This 
steady-state is achieved because ideal fluids are time- 
reversible and is thus the reverse of two equal-width jets 
incident at an angle. Due to symmetry, this reversed 
case is in turn equivalent to a jet impinging onto an infi¬ 
nite wall, a situation that surely reaches a steady-state. 
Hence, head-on collisions of ideal-fluids reach a steady- 
state where the ejecta compensate for the nonzero inci¬ 
dent momentum. 

Since this stationary steady-state impact owes its exis¬ 
tence to the fact that the flow field is entirely reversible, 
breaking reversibility by including dissipation, no matter 
how small, should destroy the stationary solution. Intu¬ 
ition would then suggest that the head-on impact of dis¬ 
sipative jets should thus produce collective motion in the 
central impact region. Numerical simulations presented 
in the current study verify this expectation by demon¬ 
strating that a wide variety of dissipational jet impacts 
produce quasi-steady-state flows whose entire flow field 
and free surface translate at a constant drift speed. This 
drift speed is shown to increase as power law relative to 
the dissipation with exponent 1.5. More surprisingly, we 
find that fluids with widely varying rheologies all collapse 
onto the same power law: the behavior of this bulk mo¬ 
tion depends only on the amount of dissipation in the 
impact, and is independent of the dissipation mechanism 
or the constitutive material of the incident jets. 

Our findings are at the intersection of two previously 
studied but disparate jet impact systems: that of directly 
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incident turbulent Newtonian streams and that of jet im¬ 
pact against fixed targets. In the case of the impact 
of turbulent jets, previous experiments identified that in 
fixed-length cavities, the impact stagnation point shifts 
toward the nozzle with lower incident momentum flux, 
whose offset depends on the incident momentum fluxes 
and nozzle separation m- These studies were focused 
on the stagnation point shifts, and reach steady state 
flows in the lab frame due to the fixed nozzle separa¬ 
tion. Our study is different because we investigate the 
offset mechanism and transient directly and are able to 
completely remove the nozzle effects. 

Regarding impacts against fixed targets, previous 
workers showed granular and high speed liquid jets as- 
toundingly both produce quite similar coherent ejecta 
sheets m- Because the flow from a high speed liquid 
impinging a fixed target is similar to a perfect fluid, this 
led to a discussion of whether or not a granular jet im¬ 
pact against a target behaves like a perfect fluid [THII]. 
Despite the similar ejecta from the target impacts, it was 
later shown that granular jet impact against a target pro¬ 
duces a large, stagnant dead-zone, and that furthermore 
inserting large dead-zone-like structures into perfect flu¬ 
ids barely affected the ejecta [12]. This dead zone, which 
is reminiscent of those found in dilute flows [IH [M] , is 
a stark contrast to the stagnation point in perfect fluid 
impacts, yet the ejecta remained generic in its presence, 
showing that the jet composition and internal structure 
have little effect on the ejecta when the impact is cold 
and incompressible. 

Directly colliding fluid streams furthermore parallel 
multiple applications in manufacturing processes. Op¬ 
posed jet micro reactors have been utilized to facilitate 
combustion and to produce highly controlled chemical 
reactions for polymer processing, nanoparticles, and the 
careful synthesis of organic compounds in the pharma¬ 
ceutical industry H H [I5HIII- Furthermore, opposing 
jets are employed in drying particles and fluid absorp¬ 
tion [181 [H] and mineral extraction [20]. The present 
findings are relevant to the mechanics by which the im¬ 
pact center shifts between nozzles in these diverse indus- 
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FIG. 1. (Color online) Two-dimensional head-on impact of two jets produces a steady-state flow with collimated ejecta and a 
steadily drifting flow configuration, (a) In two dimensions, A jet of width 2R, shown in blue, is incident against a jet of with 
2kR, shown in red. Both jets are initially at 81% packing fraction and traveling toward each other at the same speed Uo. Here, 
R = 1007?g, where Ro is the grain width, and k = 0.465. (b) Snapshot at t = 2R/Uo after the two jets directly collide. Here 
we indicate the coordinates x and y whose origin is located at the point of initial contact, (c) Snapshot at t = IQR/Uo, where 
the central saddle point, shown with a dot, is drifting toward the smaller oncoming jet with drift speed Ud- Here, Ud = 0.052[/o. 
Additionally, the impact produces 2 collimated ejecta streams with widths 2i?out and ejecta velocities with magnitude Uout 
exiting at an angle To from the horizontal. Note Tq is not quite the same as the spatial angle formed between the ejecta jet 
and the horizontal due to the fact that the ejecta streams are drifting alongside the saddle point. 


trial examples, in particular highlighting the effects of 
dissipation which have not yet been considered deeply, 
and are explored in this paper. 


II. METHODS 


In this paper, we consider the two-dimensional head- 
on collision of two dissipative jets of different sizes: one 
jet of width 2R and a smaller jet of width 2kR, where 
0 < fc < 1. The two jets are incident such that their 
centers align. Their initial velocities are both Uq in mag¬ 
nitude, and are opposite in direction so that the two 
jets collide directly (Fig. (a)). To simulate these im¬ 
pacts we used two separate numerical approaches: the 
first being a discrete particle simulation of dense gran¬ 
ular jets [n [m EH [22], and the second being the con¬ 
tinuum dynamics of colliding jets using the multi-phase, 
incompressible fluid flow solver GERRIS E51126) . 

The discrete particle simulations consist of simulat¬ 
ing the motion and collisions of hundreds of thousands 
of rigid spheres with specified coefficients of restitution 
and friction. We conduct our discrete particle simula¬ 
tions using rigid grains, in line with experiments which 
used copper and glass beads with negligible deformation 
during impact [Tunj. The discrete particle simulations 
utilize a hybrid timestep- and event-driven method de¬ 
signed to handle dense granular flow efficiently and accu¬ 
rately mull]- This method has been used in many recent 
studies of granular jet impact, and has also been quanti¬ 
tatively validated against experimental granular jet im¬ 
pact flow configurations [T^ . 

We prepare our discrete particle jets at 81% packing 


fraction, and because our simulations are in two dimen¬ 
sions, we use polydisperse grains to prevent spurious crys¬ 
tallization. The grains we use have radii uniformly dis¬ 
tributed from 0.8i?G to 1.2i?G) where Rq is the average 
grain size. In this paper, we use a grain size Rq — i?/100. 
The jet-width-ratio k, restitution and friction coefficients 
are varied and the effects on the dynamics of head-on col¬ 
lision are explored. 

The continuum simulations of colliding streams are 
conducted for liquid jet impacts, where the shear stress is 
proportional to the shear rate, for frictional fluid models 
of granular jets, where the shear stress is proportional to 
the pressure, and a shear thickening fluid whose effective 
viscosity increases with the shear rate. 

Continuum solutions in GERRIS utilize an incompress¬ 
ible volume of fluid method to solve multi-phase flow 
while providing significant flexibility for specification of 
the deviatoric stress tensor [23H25] . Ideally, we are solv¬ 
ing Cauchy’s momentum equations with incompressibil¬ 
ity 


p—=V-cr, V-u = 0 (1) 

where D /Dt = dt + u-'V is the advective derivative, p is 
the jet density, u is the fluid velocity field, and cr is the 
stress tensor. The initial conditions are 


r Uox a: < 0 

[ —Uqx X > 0 


( 2 ) 


with an initial jet free surface corresponding to the two 
impinging jets right when they make first contact. The 
boundary condition on the free surface corresponds to 
zero-stress with cr • n = 0, where n is the surface normal. 
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The boundary conditions within the jet are u = C/qX for 
X —)■ —00 and u = —Uox for x —)■ oo. 

Our stress tensor is in general given by 

(T = -pi + T (3) 

where 1 is the identity matrix, p is the pressure and t 
is the deviatoric stress. In the cases presented in this 
paper, we consider three deviatoric stress tensors: that 
of Newtonian liquid flow, where XLiq = 777 , that of a 
frictional fluid, where Tff = MP7/l7li and that of a shear 
thickening fluid tst = 7 t | 7 |^^^ 7 . Here, 7 = Vu + 
is the rate of strain tensor, and I7I = 1/7 : 7/2. 

For the Newtonian rheology, 77 is the dynamic viscosity, 
for the frictional fluid, p, is the dynamic frictional coef¬ 
ficient, and for the shear thickening fluid k is the flow 
consistency index. Speaking in particular to the fric¬ 
tional fluid, this frictional fluid model has been shown to 
reproduce simulations and experiments of dense granu¬ 
lar jet impact, and is therefore a natural choice for this 
study [nmi]. We additionally considered the Newto¬ 
nian and shear thickening fluids because their stress ten¬ 
sors are significantly different from dry granular flow. 
These differences are key in our results showing that 
widely varying materials produce similar head-on impact 
responses. 

Even though our idealized problem consists of free- 
surface jets of a single fluid type colliding in free-space, 
the volume of fluid method that GERRIS utilizes to solve 
the equations based on the boundary and initial condi¬ 
tions above requires a surrounding fluid occupying the 
free space outside the jets in our idealized problem. In 
order to minimize the effects of the surrounding fluid, we 
set its density to Psurrounding = / 5 / 25 O, and we set its devi¬ 
atoric stress to be Newtonian with viscosity ?7surrounding = 
10“^i?f7oPinter- We conducted our simulations with re¬ 
finement 26 grid points across the jet width. 

We have extensively varied the grid refinement, the 
surrounding fluid viscosity, and surrounding fluid density 
to ensure that the above parameters are in a regime where 
the effects of the surrounding fluid and grid discretization 
are minimal. Specifically, using a frictional fluid with 
p = 0.3 as a representative impact, we found that there 
was a variation in the drift speed of roughly 0.1%, and a 
variation of the energy dissipation rate of roughly 2%. 


III. RESULTS 
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FIG. 2. (Color online) Following a short transient, the central 
saddle point drifts at a relatively steady speed Ud- Normalized 
(a) saddle position and (b) corresponding drift speed versus 
time for collision between two jets at fc = 0.465 (red) and 
k — 0.25 (blue). The impact initially occurs at Xsaddie = 0. 
Following a short transient of duration ~ AR/Uo, the saddle 
point translates linearly in time. The direct measurement of 
the drift speed obtained by differentiating the saddle position. 
The drift speed starts high in the transient state and then 
fluctuates about the mean drift speed after steady drift has 
been reached. The drift speeds above are Ud/Uo = 0.052 for 
k — 0.465 and Ud/Uo ~ 0.100 for k — 0.25. 


then consider multiple continuum jet impacts and find 
that the drift behavior was extremely robust, to the ex¬ 
tent that the total energy dissipation controls the drift 
speed regardless of the significantly different dissipation 
mechanisms present in our continuum and discrete parti¬ 
cle simulations. This dependence manifests itself via the 
drift speed robustly scaling versus the dissipation rate 
with a 1.5 exponent. 


To briefly summarize our results, using discrete par¬ 
ticle simulations, we observe the resulting flow from the 
collision of two unequal-width dense jets composed of 
dry, rigid grains in two dimensions. We reveal that after 
an initial transient, the impact produces a quasi-steady 
state corresponding to a uniformly translating flow field 
and free surface at a fixed drift speed. This steady drift 
speed increases as the size of the larger jet with respect 
to the smaller jet increases, and occurs alongside momen¬ 
tum being absorbed into the drifting impact center. We 


A. Steady saddle drift 

To begin, we consider two unequal-width granular jets 
with k = 0.465, and whose grains have coefficient of resti¬ 
tution 0.9 and coefficient of friction 0.2 (see Fig. Ha)). 
Unless otherwise stated, these are the values we use in 
our representative simulation examples in our discrete 
particle jet impacts. In Fig. [I] we show early snapshots 
(a) just before impact, (b) just after impact, and (c) after 
the impact flow configuration has been well established. 
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Lab reference frame 


Comoving reference frame 




x/R x/R 

FIG. 3. (Color online) The entire flow configuration arising from the head-on collision of unequal jets drifts steadily in time, 
(a) As viewed in the fixed lab reference frame, we show horizontal velocity Ux/Uo versus position x/R along the line running 
through the centers of the jets defined by y = 0. Velocities are shown at a variety of times, ranging from very soon after 
impact (blue) to long after impact (red). The velocity fields at these different times does not coincide when viewed from the lab 
frame, (b) Horizontal velocity profile x/R along the center line as viewed in the comoving reference frame which is translating 
steadily at the drift speed Ud along with the drifting saddle. All of the profiles from (a) are included, where all of the later 
(t > 2R/Uo) snapshots directly overlap in the comoving frame. Hence, excluding the profiles obtained in the initial transient, 
the flow configuration translates at a constant drift speed Ud, remaining preserved as it translates. 


The collision produces some familiar behavior, namely 
that two collimated ejecta streams are produced as in 
the case of impact against a fixed target [3 El m HZ]. 
The collimated ejecta are traveling with velocity Uout at 
an angle 'Fg relative to the horizontal. 

However, we also find notably unfamiliar behavior in 
this system: the steady collision of the two incident 
streams produces a steadily translating mode in which 
the impact center and ejecta streams drift at a constant 
speed Ud- As a result of this drift velocity, the spatial an¬ 
gle traced by the ejecta stream relative to the horizontal 
(Fig.[T]: c)) is slightly different from Tg, the ejecta veloc¬ 
ity angle. In contrast, the more studied case of impact 
against a fixed target, there is always a fixed steady state 
where the ejecta streamlines are equal to their pathlines. 
But, in general, this is not true for head-on collisions of 
two jets, since they obtain a quasi-steady-state where the 
flow field and free surface steadily translate. 

Let us define the saddle position Xgaddie as the central 
point in the impact at the interface where the grains of 
the large jet meet the grains of the small jet. In Fig.[^c), 
this is where the blue and red jets interface along the line 
y = 0 and moving toward the smaller incident jet. The 
saddle center is initially situated at x = 0 at the point 
where the two jets first make contact (Fig.[^b)). Because 
the granular jets are dense, and therefore cold mm, 
there is negligible diffusion and mixing between the jets 
and their interface is clearly defined. In Fig. [2 (a) we 
track the saddle position over time, and find that fol¬ 
lowing a short transient, the two jets produce a steady 
central motion where the saddle position moves linearly 
in time. 

We instantaneously measure the drift speed over time 


by computing the time derivative of the saddle position 
by fitting a slope to 10 measurements of Xsaddie over a 
time frame R/Ug. These drift speed measurements are 
shown in Fig.j^b). In Fig.j^a) ,we graph only one Xsaddie 
for 10 measurements conducted in from our simulations 
so that each Xsaddie corresponds to one instantaneous 
Ud/UomFig.^h). This ensures that each Ud is indepen¬ 
dently measured from all other Ud shown. In the discrete 
particle jets, we see that there are slight fluctuations in 
the rate of motion of the saddle point about the average 
drift speed on the order of C/g/100, with an average drift 
speed Ud/Uo = 0.052 to the right when k = 0.465 and 
Ud/Uo = 0.100 when k = 0.25, where the large jet is 
pushing the small jet. Thus, smaller k produces a larger 

Ud- 

As we next show, not only does the saddle drift steadily 
in time, but the entire flow field drifts along with the 
saddle steadily in time. We measure the flow field along 
the jet centers by spatially binning grains into square 
boxes of linear dimension 4.5i?G along the center of the 
jet at y = 0 and averaging the particle speeds within each 
bin. A typical bin contains 5 particles. 

In Fig. I^a) , we show the instantaneous horizontal ve¬ 
locity field Ux along the center line defined at various 
points in time. We see the velocity in the large left jet 

^ —4i?) is Uo, and the velocity at large x in the smaller 
jet ^ 4i?) is —Uo- The velocity fields at different 
points in time as viewed from this fixed lab frame appear 
distinct from one another. 

We now enter the comoving frame whose origin is mov¬ 
ing to the right at the drift velocity Ud obtained from 
tracking the saddle point motion as in Fig. Here we 
find a dramatically simpler picture than in the lab frame. 
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Let us denote coordinates in the comoving frame with a 
bar above the variable, notation that will be continued 

through the rest of this paper, 
ordinates in this notation are 

The comoving frame co- 

X = X — Udt 

(4) 

y = y 

(5) 

Ux — Ux Ud 

(6) 


(7) 


In the comoving frame, neglecting the short transient of 
duration ^ 4i?/C/o, the velocity fields are indistinguish¬ 
able. The clear conclusion here is that when two unequal- 
width granular jets collide, they produce a quasi-steady- 
state in which the entire flow field uniformly drifts in time 
at a fixed drift speed Ud in the lab frame. In the comoving 
frame, the impact produces a steady-state, where partial 
time derivatives of the velocity field or interface locations 
(dt) vanish, and the streamlines match the pathlines, re¬ 
sulting in the spatial angle traced by the ejecta relative 
to the horizontal to match the velocity ejecta angle. 

Having established that a central drift emerges during 
impact of two head-on granular jets of unequal widths, 
we next explore the effects of jet width, and will profile 
the flow configuration as k is varied. 
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B. Characterization of post-impact flow 
configurations 

For the discrete particle jets with coefficients of resti¬ 
tution and friction of 0.9 and 0.2, we varied the jet-width- 
ratio k and observed the resulting impact configuration 
in the lab frame to provide a basic understanding of head- 
on jet impact response. This also allows us to verify that 
the simulations are producing sensible outcomes, for ex¬ 
ample, to verify that there is no central drift when the 
incident jets have equal widths. 

The impact configuration is characterized by the drift 
speed t/rf, the ejecta speed f/out, and the ejecta angle 
^ 0 - later shown, the ejecta width -Rout is determined 
uniquely from these parameters by mass conservation 
(Eq. (l5|). 

In Fig. I^a) we see how the drift speed Ud varies as k 
is varied. When the jets are the same width, there is no 
drift, as one would expect due to symmetry. When the 
small jet is 0.465 times the width of the large jet, as men¬ 
tioned earlier, there is a drift speed of Ud = 0.052{7o, and 
when the small jet is 0.125 times the width of the large 
jet, the drift speed increases to Ud = 0.17t/o- For any 
fc < 1 we found the drift emerge; there is not a thresh¬ 
old or yield k above which there is no drift and below 
which it emerges. Rather the drift speed as a function 
of k appears simply as an emergent result of unequal jet 
widths. 

In Fig.j^b) we show the dependence of the ejecta speed 
Uout on k. Here we see much less variation than we did 
for the drift speed. Evidently, at fixed jet parameters, the 


FIG. 4. (Color online) As the jets approach equal widths, 
the drift speed continuously vanishes, (a) In the lab frame, 
we show the saddle point drift speed Ud/Uo versus the ratio of 
the incident jet widths k while keeping R — IOORg and the 
restitution and friction coefficients 0.9 and 0.2, respectively 
in the discrete particle simulation (red points). We further¬ 
more include the outcome from our control volume analysis 
for comparison, where we see the analysis perform quite well 
despite its simplicity. We see that Ud decreases monotonically 
as k increases, and eventually Ud vanishes when the jets are 
equal widths at A: = 1. (b) Ejecta velocity in the lab (red) 
and comoving (blue) frames. Note that the ejecta velocity 
appears relatively constant for all k, especially in the drifting 
frame. For the grains used here, that constant is approxi¬ 
mately B = 0.73. (c) To complete the characterization of the 
flow configuration resulting from the collision we include the 
ejecta angle in both the lab and comoving frames. 


ejecta speed does not change much even though the jet 
sizes are changed dramatically. In the cases presented, 
the smaller jet varied in width by a factor of 8 while t/out 
varies by just 13% in the lab frame. 

Furthermore, we provide the ejecta angle as the fi¬ 
nal ingredient to characterize the impact flow configura¬ 
tion. We see that 'I'o decreases as k decreases, indicating 
that when the larger jet is much wider than the smaller 
one, the ejecta are more acute and are less deflected than 
the directly vertical jets obtained at the maximum deflec¬ 
tion of 90° when fc = 1. 
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C. Bulk momentum absorption 



FIG. 5. (Color online) Momentum balance incorporates 
momentum absorption into the drifting saddle. Here we show 
the momentum absorbed by the drifting saddle Psaddie as a 
fraction of the net incident momentum Bn as we vary the 
jet-width-ratio k. We see that the drifting saddle absorbs a 
considerable fraction of the net incident momentum. In the 
discrete particle jets used here, at least 15% of the net incident 
momentum is captured by the drifting saddle. 

We now show that the drifting mode absorbs a signifi¬ 
cant fraction of net incident momentum. This absorption 
corresponds to a new term in the momentum balance for 
head-on jet collisions. Let us begin by defining the mo¬ 
mentum fluxes in the three parts of the collision: Bn is 
the incident momentum flux, Pout is the momentum flux 
of the ejecta, and Psaddie is their difference: 

Bn Baut — Bgaddle (8) 

Importantly, because of the y-symmetry of our impact 
system, P in all cases above refers only to the horizontal 
momentum along the x-axis. 

In the absence of drift, such as in ideal fluid impact, 
we know that the ejecta balance the momentum equa¬ 
tion and Pout = Bn- However, in the case of granular 
jets, or drifting jets in general, we observe Psaddie enter¬ 
ing the picture: by measuring the net incident momen¬ 
tum Bn and the net momentum in the ejecta Pout in 
our simulations, we find that their difference is nonzero. 
Namely, there is momentum being deposited into the sad¬ 
dle Psaddie- 

In Fig. we show the fraction of the initial momentum 
that is absorbed by the saddle Psaddie/Bn as a function 
of the jet width ratio k. We see that indeed a significant 
fraction of the new incident momentum is being absorbed 
by the saddle, in this case roughly 15% for the range of 
k presented. 

This absorption is due to the fact that the flow field 
is uniformly translating, which in the lab frame appears 
simply as an elongation of the larger jet and a contraction 
of the smaller jet. Hence, in a unit time, the change in 
momentum due to the large jet elongation is 2pRUoUd, 
and the change due to the contraction of the smaller jet 


is 2pkRUoUd- Thus 

Biaddle = ‘2pRUoUd{l + k) (9) 

This effect is not unique to the momentum. In fact, the 
mass and energy also display this behavior, for example, 
the same approach yields a new term in mass balance due 
to the saddle motion given by Mjn — Mout = -^saddle = 
2pRUd{l-k). 


D. Drift dependence on energy dissipation 

Now we show that the main factor controlling the cen¬ 
tral drift is energy dissipation. For many jet-width-ratios, 
we vary the energy dissipation in the granular jet by 
changing the coefficient of restitution and coefficient of 
friction of the grains. We furthermore are able to explore 
the drift by varying dissipation in the impact of contin¬ 
uum fluids using the GERRIS flow solver, where frictional- 
fluid, Newtonian and shear thickening constitutive rela¬ 
tions are considered. 

The energy dissipation rate is most easily computed in 
the comoving reference frame. As mentioned, there are 
net mass, momentum, and energy transfers to the drifting 
center in quasi-steady-state found in the lab frame, and 
for simplicity we do our analysis in the comoving frame 
where a steady state is reached to simplify the analysis. 
In the comoving reference frame, we denote Bb and Bsj 
as the incident speeds of the large and small jets respec¬ 
tively, and Bout and as the ejecta speed and ejecta 
angle, respectively. These are related to their values in 
the lab frame via 


Bb = Bo - Ud 
Bs = Bo + Ud 

Bout = \J (Bout cos ^>0 - UdY + B^ot sin^ 

, -1 Bout sin ^'o 

Wq = tan *---— 

Uont COS^o - Ud 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


The jet cross sections Bout, B, and kR are the same in the 
comoving and lab reference frames. The transformations 
of the ejecta speed and angle in the lab versus the co¬ 
moving frames can be seen in for representative impacts 
in Figs.Qb) and (c). 

In the comoving reference frame, when the impact pro¬ 
duces a steady-state flow configuration, the dissipated 
energy Bdiss is equal to the energy advected in Bn sans 
the energy advected out Bout- These energy fluxes are 
given 


Bdiss — Bin Bout (1^) 

Beginning with the energy influx. Bin is simply the sum 
of the energy input from both the large and small jets 

Bin = pB(Bb + kUg) 


(15) 
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FIG. 6. (Color online) Drift speed reaches a fraction of its maximum based only on how much of the available incident 
energy is dissipated. This fraction is achieved entirely independently of all jet-width-rat ios and di ssip ation mechanisms. Shown 
in (a) is the drift speed Ud/Uo versus the normalized energy dissipation (from Eqs. (181 and (15l) for the discrete particle 
simulations and the continuum fluids. The discrete particle simulations (DPS) with coefficient of restitution (COR) 0.9 (black, 
triangles), 0.1 (dark grey, squares), and 0.05 (light gray, circles) are shown for k — 0.465. The continuum impacts are shown 
with k = 0.314 (^), k = 0.465 (-I-), and k = 0.688 (x) for the shear thickening fluid (ST, magenta), the frictional fluid (FF, 
red), and the Newtonian (Liq, blue) cases. At any given k, all of these cases produce similar drift speeds at a given energy 
dissipation rate. For comparison, we additionally include our control volume analysis (green, dashed) at A: = 0.465, which 
qualitatively captures the increasing drift speed with increasing dissipation, but is completely unable to capture other features, 
most notably the concavity, (b) Drift speed Ud normalized by the maximum drift versus the energy dissipation rate 

Ediss normalized by its maximum E^g^- Firstly, we find a dramatic collapse of the data, indicating that when the drift speed 
reaches a given fraction of its maximum dependent only on the dissipation rate and independent of the details of the highly 
varying impact materials and jet-widths-ratios. Secondly, we find that this fraction is robustly dependent on the normalized 
dissipation rate to the power 1.5. 


Next, the energy advected out by the two ejecta jets is 


^'out — 2pi?outC^out 


(16) 


We can simplify the above slightly by noting that there 
is no net mass flux, and that all of the mass entering is 
consequently leaving. Denoting the mass fluxes by M, 
we therefore have Min = -^out, giving 


R{Ub + kUs) = 2t7out -^out 


(17) 


because Mi„ = 2R{Ub + kllA) and Mout = 4pi?outC/out■ 
Throughout the above we noted that the density p of the 
jets remains constant during the impact in dense granular 
jet collisions m- _ _ 

Combining Eqs. ( |l4| ) and 0 , we obtain 

^diss = pR [C/i + kUl - (Ub + kUs)U^^t] (18) 


In Fig. we show how Ud varies based on the frac¬ 
tion of energy dissipated during impact by comparing 
the dissipated energy (18) to the input power Em- This 
dependence is displayed for 4 different cases: the discrete 
particle granular jets like those in our previous results, 
and the three continuum fluids computed in GERRIS. The 
continuum fluids we used were a frictional fluid model of 
dense granular jets mm, a continuum Newtonian im¬ 
pact, and a continuum shear thickening fluid. The con¬ 
tinuum impacts were repeated for k = 0.314, k = 0.465 


and k = 0.688 while the discrete case is only presented 
for k = 0.465. 

We varied the dissipation rate in the discrete particle 
granular jets by specifying the coefficient of restitution 
at 0.05, 0.1, and 0.9, and for each coefficient of restitu¬ 
tion, we varied the friction coefficient from 0 to 1. We 
varied the dissipation rate in the frictional fluid model 
by changing the dynamic friction coefficient p from 0.05 
to 0.7. Furthermore, we varied the dissipation rate in 
the Newtonian impact by varying the Reynolds number 
(Re) from 3 to 100 where the Reynolds number is defined 
Re = pRUo/rj. Lastly, dissipation was varied in the shear 
thickening fluid by varying the effective Reynolds num¬ 
ber Reeff = pW'I'^Uy^ jK from 6 to 100. For the frictional 
fluid and Newtonian continuum simulations, we verified 
that the dissipation rate showed reasonable behavior as 
p and rj were varied |28j . 

In Fig. [^a) we see that increasing the energy dissi¬ 
pation rate relative to the input power in the comoving 
frame, or increasing Ediss/Ei^, increases Ud- Also, in the 
asymptotic case of no dissipation, the data shows that 
the drift speed itself also vanishes entirely. Hence, in¬ 
creasing dissipation increases the drift speed, and having 
no dissipation, such as in the head-on impact of ideal 
fluids, there is no drift. Also, in Fig. [^a) we see that, 
consistent with Fig. |^a), reducing k increases Ud- 

For all of these fluids, we found that the dependence 
of Ud on Ediss/Ein all coincide and collapse together at a 
















given jet-width-ratio. The drift speed resulting from the 
head-on impact of unequal-width jets is insensitive to the 
precise type of the dissipation, rather, it only depends on 
the amount of dissipation in the jets for a given k. The 
mechanisms by which a Newtonian fluid dissipates en¬ 
ergy are different from the generalized Newtonian shear 
thickening flnid, and in turn entirely different from that 
in a discrete particle granular jet. It is remarkable how 
the system produces such similar drift speeds for such 
different jet compositions hen controlled for the amount 
of dissipation. 

E. Maximum drift speed and collapse 


that all of our head-on impacts appear to be controlled 
by these characteristic velocity and energy scales. Hence 
when the dissipation is at a given fraction of the available 
energy ifdiss/.E^“, the drift speed is at a given fraction 
of its maximum as well, independent of the 

jet-width-ratio k, and independent of the details of the 
dissipation mechanisms of the impact. Essentially, the 
emergent drift speed in the end appears to be a relatively 
simple phenomenon. The effects of the jet-width-ratio 
are easily captured based on the total available energy to 
the system to the extent where the drift speed is appar¬ 
ently determined by a self-similar equation 

Ud = C/rV(^diss/^d”T) (24) 


Here we briefly derive the maximum drift speed avail¬ 
able to the impact. Using the maximum drift speed as a 
characteristic velocity scale, we will then show that the 
scaled drift speed versus the scaled dissipation rate col¬ 
lapses across an enormous range of simulations. 

The drift speed Ud is clearly a strictly increasing func¬ 
tion in the dissipation rate, and so it is clear that the 
maximum drift speed is obtained when the dissipation is 
maximized. Here, the maximum energy dissipated 
is simply bounded by the energy supplied to the system: 


= Ein (19) 

Next, note that when drift is present, the momentum dif¬ 
ferential of the incident streams is lower in the comoving 
frame than in the lab frame. In the lab frame, the net 
incident momentum is 


Pin = 2pRU^{l-k) (20) 

whereas in the comoving frame the net incident momen¬ 
tum is 

Pin = 2pP((C7o - Udf - k{Uo + Udf) < Pin (21) 

The question at hand is how far can the dissipation push 
the comoving incident momentum differential; the an¬ 
swer is until the two jets contain equal momenta in the 
comoving frame. This corresponds to directly vertical 
ejecta streams with = 7r/2. At the maximum energy 
dissipation, where the drift speed is maximized. Pin = 0 
and 


(Po - ur^f = HUo + ur^f 


( 22 ) 


yielding a maximum drift speed 

_ I - Vfc 
Uq I -|- y/k 


(23) 


Using as the characteristic energy scale and the 
corresponding U^^^ as the characteristic velocity scale, 
we see in Fig. [^b) that our results do not only collapse 
across dissipation mechanisms for a given fc, but further 


where /(O) = 0, /(I) = I, and / is strictly increasing. 


F. Drift speed scaling relative to the dissipation 


In Fig.j^b), aside from the dramatic data collapse, we 
see that the drift speed has a power law scaling relative 
to the dissipation rate where Ud ^ E^i^. Though this 
scaling is clearly and apparently present, we do not have 
an explanation for the scaling at this time. In particular. 


in Eq. (24|. 


it appears surprisingly that Ud = P“®'’'(Pdiss/P)((Sr)^'^ 


A particularly confounding aspect of the power law 
are that straightforward analyses predict nothing like an 
exponent of 1.5. For example, in order to recover such a 
law using dimensional analysis, our analysis would have 
to yield, in some form 


r/2 


-^diss 


(25) 


However, a physically relevant dimensional quantity de¬ 
pendent on the cube of the dissipation rate is highly ob¬ 
scure, and it is unclear to us where one might arise. 

Furthermore, linear perturbations to the system, for 
example to the similarity equation Eq. (24) to infer the 


form for /, would not yield a 1.5 or 3/2 power, but rather 
a linear dependence of / on the dissipation. For all of 
these reasons, the scaling identified in Fig. |^is puzzling. 

Below we attempt to understand the results exposed 
above by presenting a very simple control volume analy¬ 
sis. As we will see, the analysis captures the broad nature 
of the effects of dissipation and unequal jet widths on the 
drift speed from our simulation, but is far from capturing 
the simple form for / as we empirically identified from 
our simulations. 


IV. CONTROL VOLUME ANALYSIS 

Here we provide a simple approximation of the drift 
speed resulting from the head on in pact of two unequal- 
width, dissipative jets as a qualitative exercise in showing 
how dissipation and an unequal jet-size-ratio can natu¬ 
rally produce a drift speed. This approximation is largely 
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independent from measurements or observations from our 
simulations, providing a from-scratch look at the effects 
of dissipation and the jet-width-ratio. 

We conduct our analysis in the drifting frame. The 
control volume analysis consists of two rough approxi¬ 
mations: one regarding the ejecta speed, and the other 
approximating the energy dissipation rate. These two 
parts are related but distinct, since the ejecta speed de¬ 
pends not only on the dissipation rate, but also the drift 
speed. 


A. Approximating dissipation rate and ejecta 
speed 


In the following analysis, we presume that most of the 
dissipation occurs in the vicinity of the smaller jet near 
the impact saddle. This corresponds to a length scale 
L ~ CkR, a velocity scale U ^ Us, and, in two dimen¬ 
sions, a volume scale oi V ~ Here, the constant C 
is assumed to be independent of the jet width ratio and 
the dissipation rate. 

The control volume analysis begins with a very rough 
approximation of the energy dissipated during impact. 
The total energy dissipation rate, shown in Eq. (18), re¬ 


lates the energy dissipation to the change in the incident 
and outgoing energy fluxes. More plainly, 


.^diss — 


^Q.dS 


(26) 


which requires 


C = 2 (31) 

The final part of the approximation concerns the ejecta 
speed. Our analysis simply approximates the ejecta 
speed in the comoving frame as 

Uout = BUo (32) 


In the above approximation, we allow B to vary when 
fluid properties vary, such as changing the Reynolds num¬ 
ber in a Newtonian impact, or the friction coefficient in a 
discrete-particle granular impact. But, we presume that 
B does not change as k is varied. This approximation 
is motivated by Fig. where we see that C/out is ap¬ 

proximately constant for the range of k tested. In that 
example with discrete particle jets whose coefficients of 
restitution and friction are respectively 0.9 and 0.2, we 
obtain B = 0.73. 

Hence, with C determined by symmetry, the analysis 
becomes dependent on the parameter B above which ac¬ 


counts for dissipation effects. Combining Eqs. (29) and 
(32), our analysis simply yields 


Ul, + kUl-{U^ + kUs)B‘^Ul = 2k{Ul-B^U^)Us (33) 

The result of the control volume analysis, shown in 
Eq. (33) is a cubic equation for the drift speed. Below we 


explicitly show the cubic equation for Ud from expanding 
Eq. (33) and recalling Eqs. (10) and ( [TI| ): 


where the contour is along a control volume surrounding 
the jet and extending far from the impact center in order 
to enclose the entire flowing region with nonzero velocity 
gradients. This control volume is fixed in space since 
the flow in the comoving frame is in steady-state. This 
becomes, noting fluid incompressibility. 


pu _ 


/v- (^u) dH=|/u.V(u2)dH (27) 

Using dimensional analysis we can approximate Ediss 
from Eq. ([27| as 


pU^-^V ^ pCkRiU^,,, - Ui)Us 


(28) 


Above we considered an additional approximation where 
A(u^) = — Ug, as this is the change in the velocity 

squared from our final velocity state relative to our dom¬ 
inant velocity scale. This, together with Eq. (18), yields 


Hi + kUi - (C7b + kUs)Ul, = Ck{Ul - Ul,)Us (29) 

For equal-width jets when k = 1, symmetry requires 
Ud = 0, and hence tJs = UsUd = 0. In this case, Eq. (29) 
reduces to 


2Ul - 2UsUl, = CUsiUi - Ul,) 


(30) 


UdV [UdVl-k 

UoJ \UoJ l + k 

+ (34) 


Because Eq. (341 has a negative discriminant for virtually 
all k and B, it in all practicality admits a unique real 
solution for Ud, and the drift speed can be determined 
relative to Uq if provided k and B. 


B. Analogy to equal dissipation in both incident 
jets 


Here we show that the energy dissipation approxima¬ 
tion in Eq. (29) is equivalent to the approximation that 


an equal amount of energy is dissipated in fluid originat¬ 
ing in the two incident jets. That is, our energy dissipa¬ 
tion approximation in Eq. [^analogous to approximating 
the energy drop across fluid originating in the large in¬ 
cident stream as equal to that in the smaller incident 
stream. 

Beginning with the small jet in the comoving frame, we 
know that the incident mass flux for particles originating 
in the small jet equals the outgoing mass flux for those 
same particles. Because the ejecta is at a speed Uout 


pkRUs — pRs,outUout 


(35) 
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where -Rs.out is the width of the ejecta stream including 
fluid originating from only the small incident jet. As 
before, we calculate the energy dissipated in the small 
jet, denoted i?s,diss) via 


.^S,diss — pkRUg pi^s^outf^out 

= pkRUsiUi - (36) 


Thus, if we set the energy dissipated in fluid from the 
large jet equal to that from the smaller jet i^B.diss = 
Es,diss, we find that this approximation yields 


.^diss — A'B^diss “t” A/B,diss — diss 


(37) 


Thus Eqs. (36) and (37) recover Eq. (29). Therefore ap¬ 


proximating the power dissipated in both jets as being 
equal is equivalent to the rough dimensional approxima¬ 
tion done earlier. 


in Fig.j^a). The analysis captures the broadest qualita¬ 
tive nature found in the simulations, namely that there is 
no drift speed in the absence of dissipation and that in¬ 
creasing the dissipation rate monotonically increases the 
drift speed. However, the analysis fails to capture the 
clear concavity found in the multiple collapsing simula¬ 
tion data sets. As such, the analysis fails to capture the 
1.5 exponent in the scaling identified in Fig. |^b). 

To its credit, the analysis does not have much tunabil- 
ity, and it nevertheless is in order the simulation results. 
However, we find that the analysis fails quantitatively. Its 
largest shortcoming is when comparing the drift speed 
versus the dissipation rate, where the analysis predicts 
the wrong concavity for Ud versus Ediss- and therefore 
does not capture the robust scaling that we identified in 
our simulations. 


V. CONCLUSION 

C. Evaluation of the control volume analysis 


We now discuss how drift emerges amid dissipation in 
our control volume analysis, as well as compare the anal¬ 
ysis to our discrete particle and continuum jet impact 
simulations. In the comparisons against previous results, 
we first compare the control volume analysis to our dis¬ 
crete particle simulations with fixed particle properties 
(or fixed B) and varied k. We then test how the analysis 
captures the effect of dissipation by fixing k and varying 
the dissipation by changing B. 

We first identify qualitative outcomes pertaining from 
the head-on collisions of two streams in the context of 
the control volume analysis from Eq. (33). Namely, we 
naturally find the presence of dissipation produces a drift 
speed. In the absence of dissipation, the jets behave as 
ideal fluids, and thus B = 1 because the ejecta speed 
must equal the incident speed. In that case, Eq. (33) 
reduces to 


Ub{Ub -U^) = kUs{Us - Ul) (38) 

whose only real solution is obtained when Ub = Us = 
Uq, and no drift is present. Decreasing B from unity 
corresponds to increasing dissipation, where Ub < Us, 
and a positive drift speed emerges. 

Second, we check how the control volume analysis per¬ 
forms when fixing B and varying k. In particular, we 
check how the analysis performs in our granular jet sim¬ 
ulations with coefficients of restitution and friction of 0.9 
and 0.2 at various k. In Fig.^a) we show the drift speed 
versus k for the granular jets as well as for the control vol¬ 
ume analysis with B = 0.73, as obtained from the granu¬ 
lar jet noting that in the comoving frame, t/out = O.TSUq 
(Fig. I^b)). We see that the analysis captures the trend 
in Ud when B is fixed and k is varied. 

The final check for our control volume analysis is when 
the jet-width-ratio is fixed at fc = 0.465 and the dissipa¬ 
tion is varied by changing B from 0 to 1. The compar¬ 
ison of the analysis to our various simulations is found 


The impact of ideal fluid jets in two dimensions is a 
classic pedagogical problem in fluid dynamics, both be¬ 
cause it is a conceptually simple generalization of the 
classical mechanics particle collision problem to fluid dy¬ 
namics and because it is analytically tractable using clas¬ 
sical conformal mapping techniques [5HH3I]. While it 
has long been appreciated that flow reversibility in ideal, 
dissipation-free jet impact dictates that the central im¬ 
pact region remains fixed in the laboratory frame, there 
has been no study focused on how dissipation affects this 
outcome. 

In this study, we investigated the effects of dissipa¬ 
tion in two dimensional head-on jet impacts, focusing on 
the head-on impact of granular jets. We identified that 
the collision produces a steadily drifting flow configura¬ 
tion which absorbs a significant fraction of the incident 
momentum into the collectively drifting impact center. 
We further showed that energy dissipation is the main 
ingredient necessary to produce the emergent drift. To 
this point, we saw that the amount of dissipation dic¬ 
tates what fraction of the maximum drift speed is reached 
across a very wide variety of fluid dissipation mechanisms 
and jet-width-ratios, and identified that the drift speed 
robustly scales relative to the energy dissipation rate with 
a 1.5 exponent. We lastly motivated a simple control vol¬ 
ume analysis of the collision which captured the general 
effect of dissipation on the drift speed, though this anal¬ 
ysis was short of explaining the scaling identified in our 
simulations. 

The findings presented above, together with ejecta in¬ 
sensitivity to internal structure in jet-target impacts be¬ 
gin to paint a picture of extreme ejecta insensitivity to 
fluid composition during dense jet impacts. In the case 
considered in this study, multiple different types of flu¬ 
ids, from discrete particle to shear thickening, all pro¬ 
duce the same drift speeds from head-on collisions when 
controlled for total dissipation; the precise form of the 
dissipation evidently does not affect the outcome of the 
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impact, rather the only determining factor is the amount 
of dissipation. Corroborating this insensitivity in dense 
jet impact in m, when impinging granular and ideal 
fluid jets against a fixed target, the internal structures 
showed very weak signatures in the ejecta. Thus, it ap¬ 
pears that the ejecta resulting from impact of dense jets 
holds a very weak signature regarding the internal struc¬ 
ture or composition of the impact. This challenges the 
traditional scattering approach in which the ejecta are 
used to infer material properties or internal structure, 
such as has been done to identify the double-helix struc¬ 
ture of DNA, the subatomic structure of Baryonic Matter 
in high energy physics, as two of several examples. 

Furthermore, when investigating the impact between 
two free streams or objects, it is natural to simplify the 
dynamics by fixing the first object and then collide the 
second against the first. Examples where this has has 
been done when studying the off-center collision of gran¬ 
ular streams m or the formation of protoplanets [31- 
133]. However, because of the presence of central drift 
in dense, dissipative impacts, the outcome of collisions 
when fixing one object is fundamentally different than 
the free-collision counterpart. This consideration can im¬ 
pact the interpretation, and perhaps even the validity, 
of some fixed-target impact experiments mimicking the 
high-speed free collision of two dense objects. For exam¬ 
ple, since the decimation of a dust aggregate is an impor¬ 
tant consideration in protoplanetary formation, the total 
force felt during impact is a crucial experimental output. 
But, in this case, fixing one object produces larger im¬ 
pact forces than free impacts because a free impact would 
reduce the collision force via the central drift. 

We would also like to highlight that the scaling behav¬ 
ior of the drift speed relative to the energy dissipation 
rate remains a mystery. Our simple approaches were not 
able to capture the exponent, and the similarity equation 
or control volume analysis presented cannot capture the 
scaling short of significant insight. In the control volume 


analysis, we assumed that the dissipation occurred over 
a length scale dependent only on k and not the energy 
dissipated. But, evidently the dissipation length scale is 
more sophisticated, and our evidence suggests that the 
relevant impact length scale varies with the dissipation 
rate. This, taken together with the fact that the de¬ 
tails of the jet constitutive material are irrelevant, lay 
the groundwork for a more involved analysis with a more 
careful accounting of the impact length scale. In particu¬ 
lar, this would manifest as a the coefficient C in Eq. (29) 
as no longer being constant, but varying based on the 
dissipation. 

Lastly, we remark that though we believe drift should 
occur in three-dimensional axisymmetric impacts, we are 
unsure if it too is insensitive to material type or results in 
a simple scaling behavior relative to the dissipation rate. 
The control volume analysis was clearly able to capture 
the qualitative nature of the central drift speed, and this 
analysis can easily be extended to show the same result in 
axisymmetric impacts. Whether or not the finding here 
that the total dissipation is significantly more important 
than the material composition is unclear, as the analysis 
was unable to entirely capture this effect, and a three- 
dimensional study with numerical diligence is required 
to determine if all of the two-dimensional findings in this 
paper persist. 
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